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We use Maxwell's equations in a sourceless, non-uniform medium with position dependent, con- 
tinuous permeability fi and permittivity e to study the wave propagation. The general form of the 
wave equation is derived and by virtue of some physical assumptions, including fi and e as functions 
of z, the equation has significantly been simplified. Finally by introducing a smooth step dielectric 
variable we solve the wave equation in the corresponding medium which is in conform with the well 
known results. Exact double-layer solution has been given in terms of the Heun functions. 



I. INTRODUCTION 



Having the electric permittivity of a system to be an isotropic and continuous spatial function has numerous 
applications in chemistry, biophysics and electronics [IHZ]. As an example, in heavy doped regions the dielectric 
constant changes with the density of impurity and so with the position. Such a region can be found in bipolar 
transistors, p — n junctions and solar cells [2J. It is also common that permittivity is assumed to be an isotropic 
spatial function that changes continuously in solvent region that allows us to formulate a computational scheme [3]- 
Other position-dependent medium can be seen in a biological membrane, like a lipid bilayer surrounded by water. 
For this system permittivity changes from a large value in the surrounding water to a lower value in the bilayer [3]. 
Also in electric double layers (EDL) due to ion accumulation and hydration in the region we face with permittivity 
variation, with effects on electric potential and interaction pressure between surfaces [5]. The main concern in such 
studies is the behaviour of a non- uniform mixed medium (a chemical solution or a p — n junction) upon an external 
static electric field using the Poisson-Boltzmann equation. This differential equation (linear or non-linear form) can 
describe the electrostatic effects extensively, ranging from a bimolecular system [H] to an electrolyte solution [7|. The 
most regular form of this equation can be written as 

V • [e (r) (r)] = -p f (r) - £ cf *?A (r) exp 

i 

In this equation e (r) is the position-dependent dielectric, i/j (r) is the electrostatic potential and Pf(r) represents 
the charge density of the medium. Further, and c°° show the charge and the concentration of ions, T is the 
temperature, Kb is the Boltzmann constant and A (r) is a factor that depends on the accessibility of a position to ions 
in the medium. 

In the present work, we study the wave propagation inside a system with the electric permittivity e (r) and the 
magnetic permeability fi (r) , for some isotropic and continuous spatial (i.e. position-dependent) functions. Without 
external sources (i.e., Pf ree — 0, J/ ree = 0) the wave equation for the electric component of the electromagnetic wave 
(EMW) propagating in the medium with e = e (r) and (j, = p, (r) is found to be 

V 2 E - ep— t = - (Vc.V) E - (E.V) Ve - V (e + //) x (V x E) (2) 
at 

in which e = lne and jl = In /i. A similar equation can be written for magnetic component of the EMW but we shall 
find B using one of the Maxwell's equation. It is obvious that with e and /i constants Eq. (2) reduces to the well 
known wave equation. A general approach toward the solution for Eq. (2) may not be possible due to the complicated 
form of the right hand side of the equation but by some simplifications we shall find exact analytical solution for 
this wave-type equation. We shall assume that the e and /i vary only in one direction which is also the direction of 
propagation. This direction throughout the paper will be z— direction. 
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Organization of the paper is as follows. In Sec. II we find the form of the wave equation in a general system of 
coordinates and through some specifications we simplify the wave equation. In Sec. Ill we introduce the smooth step 
dielectric constant and solve the wave equation accordingly. Sec. IV solves the problem of smooth double layers. Our 
conclusion is presented in Sec. V. 

II. THE WAVE EQUATION 

We start with the sourceless (pf ree = 0, J/ ree = 0) Maxwell's equations in a medium with position-dependent 
permittivity e and permeability fj, which are given by 

V.B = 0, V.D = 0, VxE = -— andVxH=— . (3) 

at at 

Herein D =eE is the displacement vector, E is the electric field, H =^B is the auxiliary magnetic field and B is the 
magnetic field. Combining the Maxwell's equations and after some manipulations one gets Eq. (2) and upon taking 
into account that e = e (z) and \i = fi (z) it becomes 



V^E - e M |J = - (e'^j E (E.V) I'k (?' + m') k x (V x E) , (4) 

in which a prime ' denotes 4-. Our further simplification is to consider that the electric and magnetic components 
of the EMW vary in z direction only. This assumption is due to the symmetry of the medium which is physically 
acceptable. Our latter equation then reduces effectively into the following 1— dimensional three equations for the 
electric components 

d 2 d 2 \ _,dE x 



d 2 d 2 \ _,dE v 



9z 2 ^dP dz ' 



/ d 2 d 2 \ 

Next we consider 

E (r,t) = E (z) e lut (8) 
in which to is the angular frequency of the wave. A substitution in (5)-(8) yields: 

/ d 2 A fi ,, n' dEi (z) 

for i = x,y. To have the third equation (i.e. Eq. (7)) satisfied we set E z = (i.e. no longitudinal component) which is 
equivalent with the propagation of the wave in z~ direction. Having symmetry with respect to x and y one can always 
rotate the system of coordinates in z direction such that E (z) aligns with one of the coordinates (say x). Hence we 
are left with only one equation in a;— direction which is expressed as 

d 2 9 \ - , „ u! dE r (z) 

^ + ^)^<*> = ^-Sr (10) 
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and the other two components of the electric field are zero i.e., E (r,i) = xE x (z) e lut . To do further investigation one 
must know the form of /i and e in terms of z. For instance, with ji = /x = cont. and e = £q = cont. one finds 



W 2 



e 2 



= (11) 



which admits 

E x = e* ikz (12) 



where k — " and the plane wave is propagating in ±z direction. Our final remark in this section is on the form of 
the wave equation: if one considers y! = reduces to the form of the standard wave equation, but owing to the form 
of e = e (z) its solution differs from the standard wave equation. 

III. SMOOTH STEP DIELECTRIC CONSTANT 

As an example let's consider fj, = K m [i — cons, and e = K e (z) e , where K e (z) stands for a smooth function of z 
given by 



AK 9 
K e (z) = K 2 — (1 - tanh (az)) 2 . 



(13) 



in which a is a positive real constant, AK = K 2 — K\, K 2 = lim^M K e (z) and Ki = lim z ^_ 00 K e (z) . The wave 
equation (9) becomes 



(^4 Ve(z)|f ' l:, = l1 



(14) 



which after defining the following new parameters 



k 2 = —^K m K\, v 2 = -^K m K 2 



and considering a new variable 



-2az 



together with a redefinition of the electric field 



E x (z) = (-0~ IV F(0: 



it turns (with ' = -4) into 



£F" + (l-2ii/) F' + 



4a 2 



(1 - 4a 2 ) v 2 y 2 -* 2 _ v 2 -k 2 



F = 0. 



Having singularities at £ — and £ = 1 suggests to replace further 

f(0 = r (e-i) p G(o 

which after some manipulation and choices 

ii/(2a- 1) 



a = 



2a 



and 



p = - ( 1 — -\/a 2 + v 2 — k 2 ) 

2 v a y 



(15) 
(16) 
(17) 

(18) 

(19) 
(20) 
(21) 
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it reduces to the following Hypergeometric differential equation (HDE) [8 



+ - \J a 2 + v 2 - k 2 | £ G' - ™ J 1 (a - \J a? + v 2 - k 2 \ 
a I 2a z \ J 



G = 0. 



a \ a 

Comparing with the standard form of the Hypergeometric DE 

i (£ - 1) G" + {(a + + 1) £ - 7] G' + aPG = 0, 

one finds 

1 



2a 



a — \J a? + v 2 — k 2 — i{y + n) 



(22) 



(23) 



(24) 



1 r 

2a 



\/ a? + v 2 — k 2 — i (y — k) 



and 



The general solution for the above HDE can be written as 

G = C 1 F(a,/3; 7 ;C) + C 2 C 1 - 7 ^(a- 7 + l,/3-7 + l;2-7;0 



(25) 



(26) 



(27) 



in which C\ and C 2 are two integration constants. Going to the original variables now, the general solution for the 
Electric field is given by 



E x (z) = C x (a,/?; 7 ; £) + 

C 2 (1 ~ 0" F (a - 7 + 1, P - 7 + 1; 2 - T, 



(28) 



where C*i = Ci (— 1) CT+P and C*2 = C 2 (— 1) <T+P+1 7 . Let us consider that the wave comes from z = —00 and goes 
toward z = +00. Also we recall that lim 2 _ i . 00 K e (z) = K 2 = cons, which implies that lim z _ i . 00 E x (z) ~ e 1 "^™^ 2 = 
e wz . Once z — »• +00 it is clear that £ = — e~ 2az — > 0, and upon knowing that F (a, f3; 7; 0) = 1, makes the limit of the 
electric field to be 



lim E x (z) = C ie wz + C 2 e~ lvz . 

z—too 

This suggests that for this choice we must set C2 — 0, which casts the solution into 

E x (z) = C x (1 - 0" F (a, 0; 7; , 



so that 



lim E x (z) = C x e lvz = E 02 e l 



(29) 

(30) 
(31) 



Herein we consider the amplitude of the transmitted wave as Eq 2 . 

The limit of z — > — oo (and consequently £ — > — oo) can be found once we apply the following properties for the 
Hypergeometric functions 



(32) 



T{a)T{j-a) v 



+ 1 - a; 



1 



Hence 



lim F(a,P; T ,0= H» ("D" C° + ^(-iff', 

2-^-00 2-^-00 1 (#) 1 (7 — m 1 fa) 1 (7 — a) 



(33) 
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which upon (28), one finds 

lim E x {z) = C 1 



-OO 

-oc 



r( 7 )r(/3-a) 
rC8)r( 7 -a) 



Now one can show that 



-P- 



(-0 



2a 
2a 



r (a) r (7 — a) 



(-0 



7> 



2a' 
''2a 



which finally yields 



lim £ s (z) = Ci 

2 — y — 00 

^ — y — 00 



r ( 7 ) r (/? - a) 
rc8)r( 7 -a)' 



r( 7 )r(a-/3) ) 
r (a) r (7 — a) 



Note that 



r( 7 )r(a-/3) 
r (a) r (7 — a) 

is the amplitude of the transmitted wave for z — > 00 and 

r( 7 )rC8- a ) 
r(/3)r( 7 -a)' 



E01 



E 'oi 



En 



En 



(34) 

(35) 
(36) 

(37) 

(38) 
(39) 



is the amplitude of the reflected wave for z — > —00. Having these definitions together with the above equations one 
can introduce the reflection and transmission coefficients of the wave as 



R 



E' m 



and 



e q1 r 08) r (a - p) 

E02 _ r (a) r (7 - a) 
E01 ~ r( 7 )r(a-/3)- 

In the limit of a sharp step dielectric i.e., a — ¥ 00 one finds 



T : 



lim 


R = 


lim 


a—yoc 




a— ^00 


lim 

a—yoo 


T = 


lim 

a— >oo 



r (a) r (/9 - a) _ k-v_ \]K m K\ - \jK m K 2 _n x -n 2 



r( 7 )r(a- j 8) 



2 k 



m + n 2 

2ni 
ni + n 2 



(40) 
(41) 

(42) 
(43) 



where n\ and 712 are the optical indices of the spaces at z < and z > 0, respectively [9]. In Fig. 1 we plot Re [E x (z)J 
in terms of 2 together with K e (z) for certain values of parameters. The smooth change in the dielectric constant and 
amplitude of the electric field are clear. 

To complete our solution we rewrite the electric field together with the magnetic field of a plane wave solution 
which propagates in positive z direction: 



E (r,i) - x£ 02 (-0" 35 (1 - 0" F (a, /?; 7; ^ 



^<, + l,« + l ;7+ l ; + (*-^ 



F(a,/3; 7 ;0 



(44) 
(45) 
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IV. SMOOTH DOUBLE LAYERS 

Another application of the general equation is to find the EMW passing through a double layer thick shell. For this 
let us consider K m = const, with 

K — K 

K e (z) = K\ H (tanh az — tanh a (z — L)) . (46) 

Here a is a constant which in the limit a — > oo, L becomes the thickness of a flat double layer dielectric of dielectric 
constant Ki located inside another medium of dielectric constant K\. Next, we rewrite the wave equation (9) in this 
matter: 



dz 



2 



k H (tanhaz - tanh a {z -L))\)E x (z) = 0, (47) 



in which k and v are defined in (15). We follow a similar change of variable given by (16) which modifies the latter 
equation as 

in which A = exp (2aL) . Now, we replace E (£) = ^ a H (£) with a — —in /2a to find 

1-te/a ^ (« 2 -^)(A-l) = 

VS7 £ 4a 2 A£(£-l)(£-i) 

which is a Heun differential equation [TO] of the form 

«/' (z) + (l + - + — ) «/ (z) + ~ 9 (z) = (50) 

\z z— 1 z — pj z(z—l){z — p) 

in which e = a + /3 — 7 — 5 + 1, with a general solution 

w (z) = CiHeunG (p, q, a, (3, 7, 5, z) + C^z 1- '* HeunG (p, q — (pS + e) (7 — 1) , j3 — 7 + 1, a — 7 + 1, 2 — 7, S, z) . (51) 

Herein C\ and Ci are two integration constants. Following the general solution we find the final solution of the wave 
equation as 

E(f ) - C, «)* HeunG fc, ^ 'ffi ' " ,0, ^, ^,<u) + (5 2 ) 

Upon considering HeunG (p,q,a, /3,7, (5,0) = 1, one easily finds that limz^-4-oo E 1 (£) = Ci exp (ikz) + C2cxp(— zkz) 

which after assuming that the wave starts from z = —00 and propagates toward z = +00 one must set C2 = and 
C\ = E03. These, therefore, imply 

E (0 = £ 03 K)^ *«mG (i, ^"^"^ .O, ^ , ^,<«) , (53) 

where -E03 is the amplitude of the electric field at the limit z — > +00. Fig. 2 displays the continuous passage of EMW 
through the double layers located at z = and z = 4. 

V. CONCLUSION 



In conclusion we have considered media with variable permeability and permittivity varying with position. The 
general form of the wave equation in such a non-uniform medium is presented. Some simplifications have been 
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considered such as, e and /x are only functions of z which is also the direction of the propagation of the possible EMW. 
In such a framework we have used the sourceless Maxwell's equations to find a wave equation for the propagation 
of the EMW. We solved the wave equation for a smooth-step, variable dielectric given the constant permeability. In 
the limits we obtain the well known reflection and transmission coefficients of the plane waves in normal incidence 
on the interface plane between two dielectrics. Our results have been presented analytically and schematically. 
Furthermore one could in principle respect not only to the position dependent permeability and permittivity but also 
a time dependent form of them which may occur in chemical solutions while some reactions are presented therein. In 
anticipation of such a case the standard wave equation would modify in the form 

<9^E $E 
V 2 E - c/i-^. = (jie + /te) — + (Jte + /xe) E (54) 

and 

„ 2 „ d 2 B dB 

V*B-e^=»e- (55) 

in which an over dot means time derivative, e = e(t) and ji = fx(t). Studying this interesting problem is our future 
plan. The applications of such extremal theories may not be so clear yet but we believe that with the fast developments 
of the new detecting methods of cancerous cells or organs in biomedical optics, such detailed theories will contribute 
to have more accurate results. Finally, in the present work we have solved the smooth double layer problem in terms 
of the Heun's functions. 
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FIG. 1: A plot of Re [E x (z)) in terms of z together with K e (z) for certain values of the parameters. The smooth changes in 
the dielectric constant and amplitude of the electric field are clear. 




L=4, co/c =1, a = 2, K m =l 
K=K= k =2, K= v =4, E„=4 



FIG. 2: The incoming EMW from z < encounters with the first layer at z — 0. Re(E x (z)) has similar structure after crossing 
the second layer at z = 4. In between, < z < 4, the oscillatory behaviour evidently changes. 



